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Abstract 

We study the neutralino sector of the Minimal Non-minimal Supersymmetric Stan- 
dard Model (MNSSM) where the fi problem of the Minimal Supersymmetric Stan- 
dard Model (MSSM) is solved without accompanying problems related with the ap- 
pearance of domain walls. In the MNSSM as in the MSSM the lightest neutralino 
can be the absolutely stable lightest supersymmetric particle (LSP) providing a good 
candidate for the cold dark matter component of the Universe. In contrast with the 
MSSM the allowed range of the mass of the lightest neutralino in the MNSSM is 
limited. We establish the theoretical upper bound on the lightest neutralino mass 
in the framework of this model and obtain an approximate solution for this mass. 



^On leave of absence from the Theory Department, ITEP, Moscow, Russia. 



1. Introduction 



The analysis of fluctuations in the cosmic microwave background (CMB) using recent 
WMAP satellite data [T] and other precise measurements [2] indicate that about 22%- 
25% of the energy density of the Universe exists in the form of stable non-baryonic, 
non-luminous matter, so called dark matter [3]. Although the microscopic composition 
of dark matter remains a mystery it is clear that it can not consist of any elementary 
particles which have been discovered so far. Thus the existence of dark matter is the 
strongest piece of evidence for physics beyond the Standard Model (SM) of electroweak 
interactions. 

The minimal supersymmetric (SUSY) standard model (MSSM) is the best motivated 
extension of the SM nowadays. Within the MSSM the quadratic divergences, which 
destabilise the scale hierarchy, are cancelled [3] and the gauge coupling unification can 
be naturally achieved [5]. If R-parity is conserved the lightest supersymmetric particle 
(LSP) in the MSSM is absolutely stable and can play the role of dark matter [6J. In most 
supersymmetric scenarios the LSP is the lightest neutralino. Since neutralinos are heavy 
weakly interacting particles they explain well the large scale structure of the Universe [7| 
and can provide the correct relic abundance of dark matter if their masses are of the order 
of the electroweak (EW) scale 

Despite these successes the MSSM suffers from the so-called /x-problem. Namely, 
the MSSM superpotential contains only one bilinear term ^.{Hd^Hu)- In order to get the 
correct pattern of electroweak symmetry breaking, the parameter /i is required to be of the 
order of the electroweak scale. While the corresponding coupling is stable under quantum 
corrections, it is rather difficult (although possible [H]) to explain within Grand Unified 
theories (GUTs) or supergravity (SUGRA) why the dimensionful parameter /i should be 
so much smaller than the Planck or Grand Unification scale. 

In the Next -to-Minimal Supersymmetric Standard Model (NMSSM) [n]^[ID], which 
contains an additional SM singlet superfield S", a symmetry forbids any bilinear terms 
in the superpotential allowing the interaction of S with the Higgs doublets Hu and Hd- 
XS{HdeHu)- At the EW scale the superfield S gets a non-zero vacuum expectation value 
((5*) = s/^/2) generating automatically an effective /x-term (/ie// = As/v^) of the re- 
quired size. There is a number of phenomenological reasons which make the NMSSM and 
its modifications quite attractive. First of all fine tuning which is needed to evade the 
LEP II Higgs mass bounds is less severe within SUSY models with an extra singlet field 
as compared with the MSSM [11]. The upper bound on the lightest Higgs boson mass in 
the singlet extensions of the MSSM was studied recently in [12J. The spectrum of Higgs 
bosons in the considered models depends on how strongly the Peccei-Quinn symmetry 
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is broken in these models [T3]-[T1]. Another nice feature is related with the electroweak 
baryogenesis which is easier to achieve in SUSY models with an extra singlet field than in 
the MSSM due to additional terms in the tree-level potential [I5]-[T6]. Recently SUSY 
models with extra singlet fields including their implications for dark matter and neutralino 
collider searches [T7| and neutrino physics [18j have been studied. 

However, the NMSSM suffers from a domain wall problem in the early Universe which 
can be avoided in the Minimal Non-minimal Supersymmetric Standard Model (MNSSM) 
as will be discussed in section 2.1. In this letter we consider the neutralino sector of 
the MNSSM. We concentrate on the mass of the lightest neutralino because it can be 
absolutely stable and therefore may play the role of the cold dark matter. We establish 
a theoretical upper bound on the lightest neutralino mass in the MNSSM which depends 
rather strongly on the parameters of the considered model. In the allowed part of the 
parameter space the mass of the lightest neutralino does not exceed 80 — 85 GeV. This 
permits to distinguish the MNSSM from the MSSM and other SUSY models with an 
extra singlet superfield at future colliders. We also find an approximate solution for the 
lightest neutralino mass. It will enable us to estimate the mass of this particle if charginos 
and Higgs bosons are discovered in the nearest future. The article is organised as follows. 
In the next section we define the MNSSM in more detail. In section 3 we examine the 
allowed range of the lightest neutralino mass in the MNSSM and in section 4 we obtain 
an approximate solution for its mass. Our results are summarised in section 5. 

2. The MNSSM 
2.1 Superpotential 

As already mentioned the NMSSM itself is not without problems. The vacuum expec- 
tation values of the Higgs fields break the Z3 symmetry in the NMSSM. This leads to 
the formation of domain walls in the early Universe [TH] which create unacceptably large 
anisotropies in the cosmic microwave background radiation [20]. In an attempt to break 
the Z3 symmetry operators suppressed by powers of the Planck scale could be intro- 
duced. But these operators give rise to a quadratically divergent tadpole contribution, 
which destabilises the mass hierarchy [21J. Dangerous operators can be eliminated if an 
invariance under or symmetries is imposed [22]-[23]. The linear term AS* in the 
superpotential which is induced in this case by high order operators is too small to upset 
the mass hierarchy but large enough to prevent the appearance of domain walls. The 
corresponding simplest extension of the MSSM is the Minimal Non-minimal Supersym- 
metric Standard Model (MNSSM) [I3], [IE], [23]-[2l]. The superpotential of the MNSSM 
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can be written as 



W; 



MNSSM 



\S{H,eH^) +^S+ WMssMif^ = 0) . 



2.2 Neutralino and chargino sectors 

The neutralino sector in SUSY models is formed by the superpartners of the neutral gauge 
and Higgs bosons. Since the sector responsible for electroweak symmetry breaking in the 
MNSSM contains an extra singlet field the neutralino sector of this model includes one 
extra component besides the four MSSM ones. This is an additional Higgsino S (singlino) 
which is the fermion component of the singlet superfield S. After the breakdown of gauge 
symmetry the Higgsino mass terms in the MNSSM Lagrangian are induced by the trilinear 
interaction XS^HdHu) in the superpotential ([1]). As a result their values are determined 
by the coupling A and the vacuum expectation values of Higgs fields. The gaugino masses 
are set by Mi and M2 which are the SU (2) and U{1)y soft gaugino mass parameters that 
break global supersymmetry. In supergravity models with uniform gaugino masses at the 
Grand Unification scale the renormalisation group flow yields a relationship between Mi 
and M2 at the EW scale, i.e. Mi ~ O.5M2 . The mixing between gauginos and Higgsinos 
is proportional to the corresponding gauge coupling and the vacuum expectation value 
of the scalar partner of the considered Higgsino. Thus after the electroweak symmetry 
breaking the superpartners of the electromagnetically neutral components of the Higgsino 
doublets and of the singlino 5* as well as the electromagnetically neutral SU (2) 
and U{1)y gauginos {W3 and B) mix forming a 5 x 5 neutralino mass matrix which in 
the interaction basis {B, W3, H^, S) reads 

/ Ml -MzSwCh MzSwSfl \ 



M^o 



\ 



Ml 

M2 

-MzSwCfs MzCwCp 

Mzswsp -Mzcwsp 





-MzSwCp MzSwSp 
MzCwCp -MzCwSp 





Xv 





71 



C/3 






Aw 




(2) 



where sw = sin 9w, cw = cos 9w, = sin/3, C/? = cos/3 and /le 
introduce tan/3 = V2/V1 and v = ^/vf + 



As 



Here we 



246 GeV, where s, Vi and V2 are the vacuum 
expectation values of S*, Hd and Hu-, respectively. 

The top-left 4x4 block of the mass matrix ([2]) contains the neutralino mass matrix 
of the MSSM where the parameter /x is replaced by /ie//- From Eq. ([2]) one can easily see 
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that the neutrahno spectrum in the MNSSM may be parametrised in terms of 



A, fieff, tan/3, Ml, M2. (3) 

The mass parameters M2 and /ig// also define the masses of the charginos, the super- 
partners of the charged gauge and Higgs bosons. Since the SM singlet superfield S is 
electromagnetically neutral it does not contribute any extra particles to the chargino 
spectrum. Consequently the chargino mass matrix and its eigenvalues remain the same 
as in the MSSM, namely 

1 



2 

m ± = - 
Xi,2 2 



eff 

(4) 



(M| + /i2 + 2M2,)2 - 4(M2/Xe// " sin 2/3y 



eff 

Unsuccessful LEP searches for SUSY particles including data collected at y/s between 
90 GeV and 209 GeV set a 95% CL lower limit on the chargino mass of about 100 GeV [23] . 
This lower bound constrains the parameter space of the MNSSM restricting the absolute 
values of the effective /i-term and M2 from below, i.e. IM2I, |/Ue//| > 90 — 100 GeV. 



3. Upper bound on the mass of lightest neutralino 

Theoretical restrictions on the masses of the neutralinos cannot be established using 
directly the neutralino mass matrix because its eigenvalues can in general be complex. In 
order to find appropriate bounds on these masses it is much more convenient to consider 
the matrix M^oMto whose eigenvalues are positive definite and equal to the absolute values 
of the neutralino masses squared. In the field basis {B, W3, H^, H^, S) the hermitian 
matrix M^oMln takes the form: 

f \Mi\^ + Mlsl^ -Mlcwsw -MzSwA* MzSwB* \ 

-Mlcwsw |M2p + M|c^ MzCwC* -MzCwD* 

-MzSwA MzCwC |/ie//P + P^ {u^-M^)Cf3Sp l^*l^effC(i ,(5) 

MzSwB -MzCwD {u"^ - Ml)cpSp |/ie//P + t^^ ^*l^effSp 

\ ^f^lffCp ""f^lff^P 



where 



_ Xv 

A = + fXeffSf3 , C = M^C/3 + f^effSfS, 

B = M*Sp + ^ieffCp, D = M*Sp + HeffCp. 
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Since the minimal eigenvalue of any hermitian matrix is less than its smallest diagonal 
element at least one neutralino in the MNSSM is always light, because the mass of the 
lightest neutralino is limited from above by the bottom-right diagonal entry of matrix 
([S]), i.e. \iTL^o\ < |z/| [2^. Therefore in contrast to the MSSM the lightest neutralino in 
the MNSSM remains light even when the SUSY breaking scale tends to infinity. 

However, the obtained theoretical bound on the lightest neutralino mass can be im- 
proved significantly. In order to get a more stringent limit on |7Ti^o| one can perform an 
unitary transformation of matrix so that M^oMIq UM^oMIqU'^ , where 



U 



(l 

1 







\ 


Cp 



C/3 





\ 
As a result we get 

-Mlcwsw 
MzSwA 

-MzSwB 




(6) 



S/3 

1 



-MlcwSw 
\M2\^ + Mlcl, 

-MzCwC 

MzCwD 




MzSwA* 
-MzCwC* 



!^ — sm 4/5 



-MzSwB* 
MzCwD* 



Ml) 



sin A(3 










where 



A 
C 



M|sin2 2/3+ |z/|2cos2 2/?^ 
Ml* sin 2p + fieff, 
M*sin 2p + fieff, 



B 
D 



M| cos2 2/?+ |i/|2sin2 2/?^ 
M* cos 2/5 , 
M* cos 2/5. 



Since we can always choose the field basis in such a way that the bottom-right 2x2 block 
of the mass matrix ([7]) becomes diagonal its two eigenvalues also restrict the mass interval 
of the lightest neutralino. In particular, the absolute value of the lightest neutralino mass 
squared has to be always less than or equal to the minimal eigenvalue /Xg of this submatrix, 
i.e. 



pI 



?r2 



-4|z/Pa2 



The value of /io decreases with increasing |/ie//|, hence reaching its maximum value, i.e. 
/Xq = min{o"^, |i^P}, when /Xe// 0. Taking the LEP bound on n^jf into account and also 



,(7) 
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the theoretical upper bound on the Yukawa couphng A which is caused by the requirement 
that the perturbation theory is vahd up to the Grand Unification scale, requiring A < 0.7, 
we find that /xg < 0.8M|, i.e. \m^o\ < 80 - 85 GeV. 

Here it is worth to notice that at large values of the effective fi-term the theoreti- 
cal restriction on |m^o| ^ tends to zero independently of the value of A. Indeed, for 
|/ie//P ^ M| we have 



l^x?r < ^ • (9) 

|/ie//|2 + 0-2 + |z/|2 



Thus in the considered limit the lightest neutralino mass is significantly smaller than Mz 
even for large values of A ~ 0.7. 



4. Approximate solution 
4.1 Characteristic equation 

The masses of the lightest neutralino can be computed numerically by solving the char- 
acteristic equation det {M^o — xJ) = 0. In the MNSSM the corresponding characteristic 
polynomial has degree 5 because the neutralino spectrum is described by a 5 x 5 mass 
matrix. After a few simple algebraic transformations we get 

det (M^o - xJ) = ( M1M2 - (Ml + M2)x + x2 ) ( x^ - (^2^^ + z/2)x+ 

+i^^f^eff sin 2/5 + M| M - X x2 + fi^^f sin 2/5x - i^M = , 



where M = Mic^+M2S^. Although one can find a numerical solution of Eq. ffTOl) for each 
set of the parameters it is worth to derive either an exact or approximate solution of 
the characteristic equation ( fTOl) to explore the dependence of the lightest neutralino mass 
on these parameters. Unfortunately, in the general case the exact solution of this equation 
is very complicated. But in the limit when one of the eigenvalues of the mass matrix ([2]) 
goes to zero one can obtain an approximate solution of Eq. (|TU|) . Indeed, if x — > we 
can ignore all higher order terms with respect to x in the characteristic equation keeping 
only the term which is proportional to x and the x-independent one. The application 
of this method is justified in the MNSSM because the mass of the lightest neutralino is 
limited from above and the upper bound on \m^o\ tends to zero with increasing |/ie//| or 
decreasing A, as argued in the previous section. Actually one can easily check that for a 
reasonable choice of the parameters (/ie//, M2 > 200 GeV, A = 0.1 — 0.7, tan P = 3—20 and 
M2 ~ 0.5Mi) the lightest neutralino mass is always significantly less than the mass of the 
second lightest one. Therefore, we can expect that the approximate solution obtained in 
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this way would describe the exact one with high accuracy in a large part of the parameter 
space. 

However, if we proceed in that way it would mean that we would allow only one 
neutralino to be light. Then the lightest neutralino mass will be consistently described 
if the four other neutralino states are considerably heavier than Mz- One can expect 
that at least three neutralino states which correspond to the superpartners of the neutral 
components of the Higgs doublets and of the neutral SU{2) gauge boson satisfy this 
requirement, because IM2I < Mz and |/ie//| < Mz are ruled out by chargino searches at 
LEP. But the mass of the neutralino state which is predominantly the superpartner of 
the U{1)y gauge boson is set by Mi which may have a value below Mz- If there are two 
light states in the neutralino spectrum then the coefficient in front of the linear term with 
respect to x in Eq. ffTUj) may be relatively small. In this case the term which is proportional 
to should be taken into account as well in order to obtain a suitable approximate 
solution for the mass of the lightest and second lightest neutralino. The inclusion of the 
quadratic term improves the agreement between the numerical and approximate solutions 
even when the second lightest neutralino is heavier than Mz- Omitting all higher order 
terms involving x" with n > 2 in the characteristic equation we find 

A>^ - Bk + C = (11) 

where 

z/^-M|/i,^^sin2/3 M| M 

^ lilff + Z/2 Ml + M2 + Ml + M2 ' ^ ^ 



M1M2 f M} M . 

B = , + „2 , .,2 - „2 , .,2 , „^ /^e//Sm2/3- 



Mi + M2 VAie// + f^lff + ^^Mi + M2 



M|i/2 



(13) 



(Mi+M2){filff + u 



z/2 / M1M2 ^ M , 



C = —o o .\ I^effsm2p-— —M'z]- (14) 

filjf + u^\Mi + M2^" M1 + M2 ^' ^ ^ 



In order to reduce the characteristic equation ( iTOl) to Eq. (jTTj) we have divided both parts 
of this equation by (/ig// + '^^)(^i + M2)- One can simplify Eq. ffTT]) even further taking 
into account that the second and last terms in Eq. ffT^ can be neglected since they are 
much smaller than unity in most of the phenomenologically allowed part of the MNSSM 
parameter space. Then the mass of the lightest neutralino can be approximated by 



Im^ol = Min <^ - 



B - VB^ - AC 



B + - AC 



(15) 
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4.2 Numerical results and discussion 

In Fig. 1 (a)-(c) we plot the numerical and approximate solutions for the lightest neu- 
tralino mass as a function of fieff, M2 and tan/5. For simplicity we assume that all param- 
eters (in]) appearing in the neutralino mass matrix are real. We also choose Mi = O.5M2 
and A = 0.7 which is the largest possible value of A that does not spoil the validity of 
perturbation theory up to the GUT scale. From Fig. 1 (a)-(b) it becomes clear that \fn^o\ 
attains its maximum at certain values of M2 and /ie//- The corresponding maximum value 
of \m^o\ is always less than the upper bound on the lightest neutralino mass derived in 
the previous section. 

As follows from Fig. 1 (a)-(c) the approximate solution (fT5l) describes the numerical 
one with relatively high accuracy even for small M2 ~ /ig// — 200 GeV. One can also see 
that the mass of the lightest neutralino may be very small or even zero for large values 
A ~ 0.7. This happens because the determinant of the neutralino mass matrix ([2]) is zero 
for a certain relation between the parameters ([3]), namely, when 

MiMa/ie// sin 2/3 = MM| . (16) 

The condition ffTUj) is fulfilled automatically when Mi ~ M2 0. It means that the 
lightest neutralino mass always vanishes when Mi and M2 go to zero. At the same time 
condition ( flGll can be satisfied at non-zero values of the soft gaugino masses. This can be 
seen in Fig. 1 (a)-(c). Since in Fig. 1 (a) and (c) we plot \Tn^o \ for non-zero values of the 
soft gaugino masses the lightest neutralino mass vanishes only once. At the same time 
in Fig. 1 (b) where we examine the dependence of \my^o\ on M2 the mass of the lightest 
neutralino vanishes twice: once for M2 = and once for a non-zero value of M2 that 
obeys Eq. ffTUj) . In the approximate solution ffT^ the vanishing of the mass of the lightest 

neutralino corresponds to the vanishing of C, which is proportional to the determinant of 

_^ det MyO 

the neutralmo mass matrix 021), i.e. C = --^ ;ttttt t-tt- 

" {filjf + u^){Mi + M2) 

Figure [2] shows the contours of the difference between the exact and the approximate 
solution ( fT5l) in the fieff-M2 parameter plane. It can be seen that this difference is smaller 
than 1 GeV in most of the phenomenologically allowed parameter space, in large regions 
even smaller than 0.1 GeV. 

Finally we would like to add that the two solutions of the reduced form of the charac- 
teristic equation (ITTl) describe with good accuracy not only the lightest neutralino mass 
but also the mass of the second lightest one if the second lightest neutralino is consider- 
ably lighter than the other states. Such a pattern of the neutralino spectrum is realised, 
for example, when Mi <C M2, /ie//- Although it is rather difficult to find any justification 
of this scenario within SUSY GUT or string inspired models it is not excluded by either 
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Figure 1: Mass of the lightest neutrahno in the MNSSM (soUd), its upper bound according to Eq. ([8]) 
(dashed) and its approximate solution according to Eq. (I15p (dotted) for A — 0.7, M\ — 0.5Af2 and (a) 
tan/3 = 5, M2 = 200 GeV, (b) tan/3 = 5, ^e// = 200 GeV, (c) ^[>'eSS = 200 GeV. 



9 



Ma [GeV] Ms [GeV] 




f^eff [GeV] /.e// [GeV] 

Figure 2: Contours of the absolute value of the difference between the mass of the lightest neutralino in 
the MNSSM and its approximate solution ((15]) in [GeV] for A = 0.7, Mi = O.5M2 and (a) tan/3 = 5, (b) 
tan/3 = 30. In the shaded region is m-± < 100 GeV. 



LEP or Tevatron searches. If z/, Mz ^ Mi in the considered hmit then the mass of the 
second hghtest neutrahno can be approximated by 

M1M2 

In this case the lightest and the second lightest neutralino are predominantly singlino and 
the superpartner of the f/(l)y gauge boson. 



(17) 



4.3 Approximate solution for decoupling limit 

With increasing effective /i-term and soft gaugino masses the lightest neutralino mass 
decreases (see Fig. 1 (a)-(b)). From Fig. 1 (a)-(c) it becomes clear that the difference 
between the numerical and approximate solutions reduces when fieff, Mi and M2 grow. 
If either /ig// or Mi and M2 are much larger than Mz, ^ C and the approximate 
solution for the lightest neutralino mass can be presented in a more simple form: 

According to Eq.( fT8i) the mass of the lightest neutralino is inversely proportional to the 
effective //-term. It vanishes when A tends to zero. In the limit A ^ the equations 
for the extrema of the Higgs boson effective potential which determines the position of 
the physical vacuum imply that the vacuum expectation value of the singlet field rises as 
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Mz/\- In other words the correct breakdown of electroweak symmetry breaking requires 
/ie// to remain constant when A goes to zero. As a result from Eq. (|T8l) it follows that 
the mass of the lightest neutralino is proportional to at small values of A. At this 
point the approximate solution ffTSl) improves the theoretical restriction on the lightest 
neutralino mass derived in the previous section. This is because at small values of A the 
upper bound ([9]) is proportional to A. 

From Eq. (fT8|) one can also see that the mass of the lightest neutralino decreases 
when tan/? grows. The numerical results of our analysis summarised in Fig. 1 (a)- 
(c) confirm that becomes smaller when tan/3 raises from 3 to 10. However, if 

tan/3 > C = ^ — -^7-^) Eq- dlHl) does not provide an appropriate description of the 

lightest neutralino mass. Indeed, in accordance with Eq. ( ITHl) the mass of the light- 
est neutralino vanishes at large values of tan/? while Fig. 1 (c) demonstrates that \'m^o\ 
approaches to some constant non-zero value with raising of tan (3. More accurate consid- 
eration of the approximate solution ffT^ allows to reproduce the asymptotic behaviour of 
the lightest neutralino mass at fj,eff,M2,Mi 3> Mz and at very large tan/? ^ C- It is 
given by 



M 



M1M2 



(19) 



So once again the approximate solution ( iTSll improves the theoretical restriction on the 
lightest neutralino mass because the upper limit ([8])-([9]) on \m^o \ obtained before depends 
rather weakly on tan (3. 



5. Conclusions 

In this letter we have examined the theoretical restrictions on the lightest neutralino mass 
within the Minimal Non-minimal Supersymmetric Standard Model. In order to derive 
the appropriate upper bound we consider the hermitian matrix M^oMIq where M^o is 
the neutralino mass matrix. The eigenvalues of this matrix are the absolute values of 
the neutralino masses squared. Therefore all eigenvalues of M^oMIq are positive definite. 
Using the theorem that the smallest diagonal element of a hermitian matrix is always 
larger than the minimal eigenvalue of this matrix we establish an upper bound on the mass 
of the lightest neutralino m^o in the MNSSM. The direct application of this theorem leads 
to the conclusion that \m^o\ has to be always less than \X\v/V2. A more stringent limit 
on the lightest neutralino mass can be obtained by applying an unitary transformation to 
the matrix McoML. As a result we have found that Im^ol does not exceed 80 — 85 GeV. 
The corresponding upper bound depends rather strongly on the effective /i-term |/ie//| 
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which is generated after the electroweak symmetry breaking. At large values of |/Ue//| 
the upper limit on |?77.^o| goes to zero so that the mass interval of the lightest neutralino 
shrinks drastically. 

Assuming that \m^o\ is considerably less than the masses of the other neutralino 
states we have derived an approximate solution for the lightest neutralino mass. The 
obtained solution describes the numerical one with high accuracy in a large region of the 
phenomenologically allowed part of the MNSSM parameter space. Our numerical analysis 
and analytic considerations show that m^o decreases with increasing tan (3 and decreasing 
coupling A. At small values of A the mass of the lightest neutralino is proportional to 
A^. The lightest neutralino mass also decreases with increasing /ie//, Mi, and M2. We 
have argued that at large values of the effective /i-term m^o is inversely proportional to 
fj'eff- In the allowed part of the parameter space the lightest neutralino is predominantly 
singlino that makes its direct observation at future colliders challenging. In forthcoming 
publications we plan to consider the potential discovery of such a neutralino at the LHC 
and ILC. 

In summary, the obtained theoretical restriction on the lightest neutralino mass allows 
to discriminate the MNSSM from other SUSY models where the mass of the lightest 
neutralino is not limited from above. If no light neutralino is detected at future colliders 
the MNSSM will be ruled out. 
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